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We  proposed  in  Phase  I  of  the  current  work  to  compuationally  model 
a  capillary  discharge  x-ray  laser  driver  which  will  use  Ne-like 
argon  plasma  as  the  active  medium  and  which  will  explore 
collisional  excitation  to  pump  upper  lasing  levels.  During  Phase  I 
we  committed  to  perform  atomic  physics  calculations  regarding 
possible  lasing  level  pairs  in  Ne-like  Ar,  then  kinetic 
calculations  to  find  level  populations  and  estimate  a  gain.  We  also 
expected  to  elaborate  MHD  simulation  of  plasma  behavior  in  a 
capillary  discharge  and  test  the  code  by  fitting  it  with  known 
capillary  discharge  experiments. 

In  Phase  II  of  the  project  we  will  build  the  capillary  discharge  x- 
ray  laser  driver  according  to  the  design  of  the  Phase  I.  We  will 
test  the  discharge  for  lasing  m  proposed  lines. 

Phase  I  was  very  successful.  We  have  performed  more  than  we 
expected  in  this  Phase  of  the  current  Proposal. 


I.  ATOMIC  PHYSICS  CALCULATIONS. 

First,  making  a  literature  survey  we  found  that  for  3-3  transitions 
of  Ne-like  Ar  there  are  still  no  sufficiently  reliable  atomic 
calculations  for  the  radiative  transition  probabilities.  We  needed 
values  of  radiative  decay  rates  to  find  out  which  energy  level 
pairs  are  candidates  for  lasing  and  eventually  for  a  gain 
calculation.  It  was  found  that  some  high-lying  level  positions  have 
also  been  calculated  with  marginal  accuracy  to  identify  transitions 
between  them  in  a  real  plasma. 

We  decided  to  perform  the  required  atomic  physics  calculation 
ourselves.  These  calculations  include  a  determination  of  the  energy  !  i-J 

level  positions  of  a  given  multicharged  ions  as  well  as  a  _ ^ 

calculation  of  the  radiative  transition  rates  and  the  branch  * 

1  I  .  .1:0 rf 


rLCT, J&D 


A.  iiljb.i  tv 


Aval  j  ■  1 
V'-C-  i 


u 


I 


coefficients  for  the  strongest  and  most  important  transitions  in  a 
multicharged  ion. 

We  performed  these  calculations  for  Ne-like  Ar,  Hi-like  Kr  and  Ni- 
like  Mo  ions.  For  Ne-like  Ar  all  36  levels  with  n=3  were  included 
in  the  calculations,  and  for  Ni-like  ions  we  are  able  to  include 
about  100  general  levels  of  the  n=4  configuration. 

The  calculations  were  performed  by  the  Relativistic  Perturbation 
Theory  Method  with  the  Model  Potential  of  zero  approximation 
( RPTMP )  .  This  method  is  well  developed.  It  has  been  shown  to  be 
reliable  for  calculating  energy  levels  for  various  isoelectronic 
sequences.  RPTMP  is  especially  successful  in  the  calculation  of 
rate  constants  of  processes  involving  multicharged  ions,  photons, 
and  electrons.  A  detailed  description  of  method  is  given  in  the 
Mid-term  Technical  Report  of  current  project.  Here  we  outline  the 
foundations  of  our  approach  which  is  called  the  energy  approach 
[1,2]. 

In  a  nonrelativistic  theory  a  field  form  procedure  is  connected 
with  a  diagonalization  of  the  secular  matrix.  It  is  a  standard 
method  for  a  calculation  of  an  energy  shift  AE  of  degenerate  atomic 
states.  Our  approach  is  based  on  the  Gell-Mann  and  Low  adiabatic 
formula  and  on  a  perturbation  expansion  of  a  scattering  matrix.  We 
use  an  analogous  scheme  with  an  electroth  electrodynamics  scattering  matrix. 
The  perturbation  expansion  for  the  energy  has  only  even  powers.  The 
contributions  are  usually  presented  Jn  the  form  of  Feynman 
diagrams.  An  important  advantage  of  the  approach  is  the  effective 
inclusion  of  high-order  perturbation  theory  (PT)  corrections. 

The  method  was  tested  to  compare  calculations  with  experimental 
measurements  for  transitions  in  Ne-like  ions  and  was  found  to  have 
accuracy  better  than  10  30  cm'1  for  level  positions  for  elements  with 
Z  <  50. 

A  more  detailed  description  of  the  calculation  scheme  as  well  as 
the  results  of  our  calculations  in  the  form  of  tables  for  three 
mentioned  ions  was  presented  in  the  mid-term  report  on  current 
proposal . 

Here  we  discuss  only  the  transitions  in  Ne-like  Ar,  and  Ni-like 
ions  of  Kr  and  Mo  we  found  to  be  most  promising  for  lasing  in  a 
capillary  discharge. 


Ne-like  Ar. 

In  Ne-like  Ar  these  transitions  are: 

2 P1/2  3p3/2  [J=2]  —  2pi/2  3Sj [J=l],  X  =69.76  nm, 
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2P3/2 

3P3/2 

[J=0] 

“  2pw 

3s,a 

X  =43.07 

nm. 

2Ptn 

3Pj/2 

[J=0] 

"  2p„2 

3Sl/2 

X  =71.25 

nm, 

2Pt/2 

3  Pm 

[J=l] 

_  2p1/2 

3sm 

X  =72.75 

nm. 

Special  attention  should  also  be  paid  to  the  level  2s2p43s  (J=0) 
which  corresponds  to  the  excitation  of  the  2s  electron.  The  energy 
of  this  level  excitation  in  Ne-like  Ar  is  comparatively  high  (330 
eV) ;  nevertheless,  the  electron  collision  rate  for  its  population 
is  high  enough  to  provide  an  efficient  inversion  population 
relative  to  the  low  active  level  2pV23sI/2( J=1 ) .  So  the  search  for 
lasing  at  the  transition 

2S tn  3sm  [J=0]  -  2p3/2  3s, ,2  X  =15.73  nm, 
in  Ne-like  Ar  will  be  of  significant  interest. 


Ni-like  ions. 

A  Ni-like  ion  is  much  more  complex  to  compute  than  Ne-like  one. 
Nevertheless,  the  RPTMP  method  was  also  successfully  applied  ►<-< 
such  a  complex  ion. 


Ni-like  Kr. 

The  numerical  results  of  calculation  for  states  3p6  3d9  4s,  4p,  4d, 
4f  and  3p5  3d10  4s,  4p,  4d,  4f  of  Ni-like  krypton  are  given  in  Table 
2  of  the  mid-term  technical  report  for  Phase  I  of  the  current 
Project.  In  this  table  we  calculated  position,  lifetime,  and  a 
probability  of  radiative  decay  to  the  ground  state  for  each  of  the 
92  levels  of  Ni-like  Kr.  The  strongest  radiative  transition  channel 
as  well  as  the  wavelength  of  this  transition  and  a  branch  factor 
were  also  calculated  for  each  of  these  92  levels. 

For  Ni-like  ion,  levels  3dw  4p1/2  [J=l],  3dw  4p,/2  [J=l],  and  3d„2 
4 p3/2  [J=l]  can  serve  as  low  active  levels. 

Four  laser  transitions  were  found  as  most  possible  in  Ni-like  Kr: 


3^5  n 

^^5/2 

[J=0] 

-  3d, ,2 

^Pl/2 

[J=l], 

X  =31.21 

nm. 

3d3n 

4d3/2 

[J=0] 

3d, /2 

4Pi/2 

X  =42.56 

nm. 

3^3  n 

4d3/2 

[J=l] 

"  3dw 

*Pm 

X  =41.50 

nm, 

3^3/2 

4d5/2 

[J=2] 

~  3d, ,2 

^  P3/2 

X  =41.78 

nm . 

The  above  outlined  approach  gives  us  a  quick  method  for  evaluating 
laser  possibilities  for  any  ion.  Nevertheless,  it  is  necessary  to 
look  completely  through  the  transition  probability  matrix  to 
determine  the  second  and  third  most  important  decay  channels.  Thus 
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we  are  able  to  understand  that  in  Kr  IX  the  decay  from  the  high 
active  levels  is  almost  equally  distributed  between  three  low 
active  levels.  The  most  favorable  cases  are  for  3d,7  4dJ/2  [J=0]  and 
for  3d5/2  4d5/2  [J=0]  high  active  states.  They  have  a  relatively  large 
branch  coefficient  for  the  transitions  to  low  active  level.  A  high 
active  level  3p3/2  4p3/2  (J=0]  is  also  shown  which  corresponds  to  fast 
3p  subshell  excitation. 

It  was  found  that  at  the  same  electron  temperature,  the  collision 
excitation  rate  of  this  level  is  approximately  the  same  as  the 
excitation  rate  of  both  of  the  above-mentioned  levels  with  J=0. 
However,  there  is  no  well-separated  dominant  channel  for  a 
radiative  transition  into  the  fast  releasing  low  active  level. 
Nevertheless,  this  level  may  influence  populations  of  upper  laser 
levels  in  the  plasma. 


Ni-like  Mo  ion. 

Calculations  similar  to  those  for  Kr  were  performed  for  Ni-like  Mo. 
The  Mo  plasma  could  serve  as  an  active  medium  in  evacuated 
capillary  discharge  with  Mo  electrodes  or  a  capillary  filled  with 
a  gaseous  Mo  compound. 

The  following  transitions  were  found  to  be  promising  for  lasing  in 
Ni-like  Mo. 


‘5/2 

4  ^5/2 

[  J=0  ]  -  3d37 

4Pl/2 

X  =17.98 

nm. 

*3/2 

4  <*5/2 

[J=2]  -  3d,, 

4P3/2 

X  =25.15 

nm. 

‘3/2 

4  <*3/2 

[J=l]  -  3d37 

4Pl/2 

[J=l], 

X  =25.37 

nm. 

The  first  two  should  have  approximately  the  same  order  of 
intensities,  and  the  last  one  is  less  intense.  In  Ni-like  Mo  two 
very  high  positioned  levels  3p37  4p3/2  [J=0]  and  3p1/2  4p,,7  [J=0]  are 
also  strongly  excited  by  electron  collisions  and  are  promising  to 
serve  as  upper  laser  levels.  According  to  our  calculations  they  are 
well  populated  by  electron  collisions  and  have  rather  high  branch 
factors  for  transitions  to  lower  active  levels.  Thus,  the  short- 
wavelength  lines 

Pi/2  4Pi/2  [ j~°]  -  3d3/2  4Pi n  ^  =6.80  nm, 

3p„2  4 P.i/2  [J=0]  ~  3d3/2  4p1/2  [J=l]  ,  X  =7.35  nm , 

could  be  of  interest  for  lasing  and  should  be  experimentally 
investigated  in  more  detail  by  spectroscopic  methods. 

The  level  3d3/2  4d3/2  [J=l]  is  not  included  in  this  scheme  because  it 
was  found  to  have  too  small  branch  coefficients  for  transitions  to 
the  low  active  levels. 
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II.  PLASMA  KINETICS  CALCULATIONS. 


Calculation  of  the  kinetics  of  level  excitation  and  deexcitation 
was  performed  to  find  level  populations  and  a  gain  in  proposed 
lines.  For  Ne-like  ions  this  resulted  in  the  creation  of  a  computer 
code  named  "POPULAT. ICN" .  The  code  calculates  (for  a  given  electron 
density  and  temperature)  the  populations  of  levels  involved  into 
consideration  versus  time.  General  rate  equations  in  the  following 
form  were  used  to  calculate  the  population  of  specific  level  n : 

dN 

+N/N1  {Rir*Rcr  +  Rde)  1  on 


Here  N,  =  electron  density,  cm3;  Nk  and  Nn  -  density  of  population 
of  levels  n  and  k,  cm'3;  =  radiative  decay  rates,  sec  x;Ckn,  C = 
electron  collisional  deexcitation  rates  from  an  arbitrary  level  k 
to  lower  specific  level  n  and  vice  versa,  cm3sec‘.  The  escape 
factor  G(t)  (where  r  is  optical  depth  of  plasma)  is  introduced  here 
to  account  for  the  effect  of  reabsorption  [3,4].  I  and  R  are 
ionization  and  recombination  rate  coefficients  for  excitation  from 
adjacent  species  (Na-like  and  F-like  ions,  respectively),  cm’sec 1 . 
The  subscripts  rr,  cr ,  de  and  stand  for  radiative,  collisional  and 
dielectronic  recombination  rates. 

At  present  we  have  completed  the  calculation  of  level  populations 
for  Ne-like  argon.  37  levels  of  Ne-like  Ar  are  included  in  the 
calculations,  including  the  ground  state,  all  26  excited  states  of 
the  2p531  (l=s,p,d)  configuration,  and  all  10  states  of  2s2p63l 
(l=s,p,d)  configuration.  Collisional  excitation  and  deexcitation 
rates  for  Ne-like  Ar  were  also  calculated  using  the  above  described 
RPTMP  method.  This  was  very  time  consuming  but  an  important  part  of 
the  current  work.  The  rates  of  excitation  and  deexcitation  (in  a 
plasma  with  parameters  of  interest  for  lasing  in  proposed  lines) 
turned  out  to  be  much  higher  than  those  between  different  ion 
species.  Taking  into  account  a  long  lifetime  for  a  capillary 
discharge  plasma  (compared  with  the  ionization  balance  change 
time),  we  facilitated  the  first  run  of  kinetics  calculations  by 
omitting  the  influence  of  neighboring  species  on  a  population. 
Nevertheless,  the  inner-shell  ionization  of  Na-like  ions  (the 
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process  creating  excited  states  of  Ne-like  Ar)  was  found  to  compete 
with  collisional  excitation  from  the  Ne-like  ground  state  and  will 
be  included  in  a  more  detailed  kinetics  calculation  scheme,  which 
we  plan  to  run  in  the  Phase  II. 


The  present  code  calculates  temporal  progress  of  population,  with 
logarithmic  time  increments.  We  assumed  for  initial  conditions  {at 
t=0)  that  a  population  of  the  ground  state  of  Ne-like  Ar  is  unity 
(100%)  and  all  other  states  are  empty  (Nk=0,  k^O)  .  This  assumption 
may  not  be  valid  in  a  real  plasma  but  it  does  not  influence  the 
final  result,  which  is  the  same  for  different  initial  population 
distributions. 


In  Table  1  we  represent  the  calculated  level  populations  for  Ne- 
like  Ar  in  a  plasma  for  electron  densities  Nt  =  1018cm'3  and  lO’Vm1, 
and  temperatures  Te  =  30  eV  and  60  eV.  At  these  temperatures  an 
abundance  of  Ne-like  ion  in  a  steady-state  plasma  with  an 
established  ionization  equilibrium  was  found  to  be  10%  and  80% 
respectively.  According  to  our  calculations,  the  time  required  to 
establish  populations  within  Ne-like  ions  for  these  plasma 
parameters  is  much  less  than  the  time  for  establishment  of  an 
ionisation  equilibrium.  Generally  this  time  is  found  to  be  less 
than  several  nsec.  After  this  time  a  permanent  inversion  population 
between  several  3p  and  3s  levels  is  established.  As  an  example,  the 
temporal  history  of  populations  is  shown  in  Fig.  10,  11  for  a 
temperature  of  60  eV  and  densities  of  1018  cm3,  and  101S  cm3. 
Numeration  of  levels  in  Table  1  and  in  Figs  1,2  correspond  to  those 
in  Table  1  of  the  Phase  I  Mid-term  Report.  Levels  2  and  4  represent 
lower  laser  2p3s  active  levels.  Levels  10,  11,  12,  14  and  28 
represent  upper  2p3p  levels.  As  was  expected,  in  a  plasma  with 
higher  electron  density  the  populations  are  established  faster. 


As  an  example,  time  history  of  population  of  these  7  levels  is 
presented  in  Fig.  1  and  2  for  two  different  electron  densities  in 
a  capillary  discharge  plasma. 


Final  populations  for  different  electron  densities  and  temperatures 
are  shown  in  Table  1.  In  this  table  we  keep  the  same  level 
designations  as  in  Table  1  of  the  mid-term  report. 
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Table  1.  Populations  N„  of  levels  in  quasi-steady-state  Ne-like  Ar 
plasma.  Populations  presented  as  a  fraction  of  total  population  of 
Ne-like  Ar  in  plasma  for  different  electron  densities  Nc(in  1018cm3) 
and  temperatures  Te  (in  eV)  .  Level  0  represents  a  ground  state, 
levels  2  and  4  represent  lower  possible  active  levels,  and  levels 
10,  11,  15  and  28  represent  upper  possible  laser  active  levels. 


Tc 

(eV) 

30 

60 

30 

60 

K 

(1018  cm3) 

1 

1 

10 

10 

§ 

Level 

J 

0 

2sJ2p6 

0 

.9965 

.6874 

.9969 

.74 

37 

2 

2p3/23s1/2 

1 

.  17 

E-05 

.  13 

E-0  3 

.97 

E-05 

.67 

E-03 

4 

2Pi/23sin 

1 

.  57 

E-06 

.40 

E-04 

.  35 

E-05 

.23 

E-03 

10 

2p,/23pl/2 

1 

.44 

E-05 

.30 

E-03 

.36 

E-04 

.23 

E-02 

11 

2Pi/23Pi/? 

0 

.  31 

E-05 

.  13 

E-03 

.25 

E-04 

.  13 

E-02 

12 

2Pi/2  3p3/2 

2 

.  82 

E-05 

.  56 

E-03 

.56 

E-04 

.  37 

E-02 

15 

2p3/23P3/2 

0 

.29 

E-04 

.24 

E-02 

.71 

E-04 

.  54 

E-02 

28 

2sin3Sin 

0 

.53 

E-06 

.59 

E-04 

.41 

E-05 

.53 

E-03 

According  to  the  kinetics  calculations,  the  most  favorable 
conditions  for  inversion  population  between  upper  levels  10,  11,  12 
and  lower  level  4  occur  at  electron  temperatures  Te  =  40-50  eV  and 
densities  around  10,9cm'3.  For  the  shorter  wavelength  transition  28-2 
the  temperature  should  be  slightly  higher  (60+  eV)  and  the 
inversion  is  generally  less. 

This  code  is  currently  being  modified  to  account  for  an  excitation 
from  neighboring  species,  mainly  from  Na-like  ion  for  our 
temperature  range  (both  from  ground  and  excited  states) .  It  seems 
to  be  a  process  contributing  essentially  to  both  upper  and  lower 
active  levels  of  Ne-like  ions  at  lower  electron  temperatures 
(around  30-40  eV) . 

Kinetics  calculation  for  Ni-like  ions  turned  out  to  be  much  more 
complicated  than  for  Ne-like  Ar  mainly  due  to  a  limited  computer 
power.  We  did  not  plan  to  perform  calculations  for  Ni-like  ions  in 
Phase  I.  Howerver,  we  did  get  the  work  started.  We  expect  to 
complete  these  calculations  in  Phase  II,  as  Ni-like  Ar  is  very 
attractive  lasder  medium,  requiring  lower  temperatures  to  populate 
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upper  laser  levels  yet  having  shorter  wavelength  transitions. 


For  the  calculated  populations  in  Ne-like  Ar  we  also  estima.  d  a 
small  signal  gain  and  possible  amplification  for  different  lines  in 
the  proposed  geometry  (2x100  mm)  of  a  capillary  discharge  x-ray 
laser  driver  (see  Table  2.) 


Table  2.  Estimated  gain  (in  cm1)  for  different  transitions  in  Ne- 
like  Ar  in  the  2x100  mm  capillary  discharge  plasma,  and  related 
amplification  of  spontaneous  emission. 


Te  (eV) 
Ne  (101* 

cm'3) 

30 

1 

60 

1 

30 

10 

60 

10 

# 

Trans i 
tion 

Wavelength 

(nm) 

Gain  G  (cm1)  ,  GL 
related  amplification 

product  and 
I/I0  =  exp  (GL) 

1 

12-4 

69.76 

G 

GL 

exp(GL) 

. 57E-3 
. 57E-2 
1.006 

.  58E-1 
.  58E  0 
1.786 

.  79E-2 
.  79E-1 
1.082 

. 12E+1 
. 12E+2 
1.6E5 

2 

11-4 

71.25 

G 

GL 

exp(GL) 

. 12E-3 
. 12E-2 
1.001 

.  68E-2 
. 68E- 1 
1.070 

. 2  IE-2 
. 21E-1 
1.021 

.  25E  0 
.  25E+1 
12 . 18 

3 

10-4 

72.75 

G 

GL 

exp (GL) 

. 15E-3 
. 15E-2 
1.002 

. 19E-1 
.  19E  0 
1.209 

. 2  5E-2 
. 25E-1 

1 . 025 

.  39E  0 
. 39E+1 
49.4 

4 

14-2 

43.07 

G 

GL 

exp (GL) 

. 51E-3 
. 51E-2 
1.005 

.  40E-1 
.  40E  0 
1.492 

. 28E-2 
.  28E-1 
1.028 

.  48E  0 
.  48E+1 
121.5 

5 

28-2 

15.8 

G 

GL 

exp (CL) 

35E-6 
-. 35E-5 

Q  9  Q  Q  q 

. 4  IE-3 
. 41E-2 
1.004 

. 17E-4 
. 17E-3 

1 . 002 

. 13E-1 
.  13E  0 

1 .  14 

Calculations  of  gain  for  other  temperatures  and  densities  were  also 
performed  (not  shown  in  the  Table  2).  It  was  found  that  a  gain  for 
the  strongest  transitions  12-4  and  10-4  reachs  a  maximum  at 
electron  temperatures  around  50-60  eV  and  densities  lO19  cm1.  At 
higher  densities  gain  decreases  essentially  as  a  consequence  of 
collisional  mixing  of  upper  and  lower  active  levels  as  well  as  due 
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to  a  reabsorption  of  a  radiation  from  the  lower  active  levels  in  a 
capillary  plasma.  The  last  effect  plays  major  role  in  controlling 
gain  in  our  capillary  discharge  plasma  at  high  densities,  as  plasma 
becomes  optically  thick  in  Ne-like  Ar  resonance  lines. 


As  can  be  seen  from  Table  2,  the  GL  product  for  the  proposed  x-ray 
laser  driver  could  be  high  enough  for  10-cm  plasma  column.  For  the 
strongest  transition  12-4  GL  reachs  the  value  12,  which  is  usually 
considered  close  to  a  saturation  regime.  The  next  strongest 
transition  is  the  short-wavelength  14-2  transition,  which  shows 
possible  a  GL  value  4.8  and  amplification  more  than  IQ2.  The  two 
transitions  10-4  and  11-4  show  moderate  GL=3.9  and  2.5, 
respectively.  This  results  in  moderate  amplification  of  50  and  12, 
which,  nevertheless  is  still  possible  to  observe  experimentally.  At 
last,  the  shortest-wavelength  transition,  28-2  (15.8  nm)  has  a 
negligible  amplification. 


For  the  Ar  plasma  with  an  ion  temperature  of  50-60  eV  the  Doppler 
width  lies  in  a  range  AX/X  =  10J.  A  low-  and  medium-resolution 
spectrometer  will  measure  intensity  integrated  over  an  instrumental 
profile.  To  measure  the  gain  properly,  an  integral  equation  for  a 
given  instrumental  factor  and  line  shape  should  be  solved. 
Fortunately,  it  was  found  [5]  that  for  any  sharply  peaked  line 
profile,  including  Doppler,  the  following  analytical  dependance  of 
measured  intensity  on  the  gain  serves  as  a  good  approximation  of 
the  exact  solution  for  this  equation: 


T_j  (exp  (GL)  -1)  3/" 
G3/2  (L  exp  {GL)  )  1/2 


Here  I  =  measured  line  intensity,  G  =  gain,  L  =  length  of 
amplifying  media,  J  =  intensity  of  plasma  spontaneous  emission  per 
unit  length.  This  dependance  was  found  to  be  the  most  correct  and 
widely  used  for  ar.  experimental  gain  measurement.  We  will  use  this 
equation  to  derive  a  gain  in  our  experiment. 


3.  CAPILLARY  DISCHARGE  PLASMA  MHD  SIMULATION 

To  determine  the  plasma  temperature,  density,  velocity  and  ion 
abundances,  a  magnetohydrodynamic  (MHD)  code  (simulating  capillary 
discharge  plasma)  was  created.  Actually  we  started  with  an 
implementation  of  the  "PINCH"  MHD  code  developed  at  Los-Ala.t.oo 
National  Laboratory  for  high-curren  deuterium  Z-pinch  plasma 
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modelling.  We  quickly  discovered  that  the  code  does  not  describe  a 
capillary  discharge  plasma  at  all.  We  tried  to  simulate  the 
available  capillary  discharge  experiments  but  we  could  not  get 
agreement  between  temperature  and  density  measured  in  these 
experiments  and  values  predicted  by  "PINCH”.  This  was  mainly  due  to 
the  fact  that  Lue  capillary  wall  processes  were  not  accounted  for 
in  the  code.  Modifications  of  the  code  did  not  help  much,  so  we 
wrote  our  own  MHD  code  for  a  capillary  discharge  plasma  geometry. 

The  one-dimensional  .capillary  discharge  plasma  radius  dependent) 
MHD  code  was  crea*'~J.  The  code  includes  a  standard  set  of  plasma 
MHD  equations  a*-  co  mass,  momentum  and  energy  conservation  (where 
energy  is  the  sum  of  thermal,  kinetic  and  magnetic  energy),  and 
Faraday's  and  Ohm's  law  equations.  These  equations  are  coupled  with 
radiation  energy  loss  and  electron  thermal  conductivity  loss 
equations . 

A  crank-Nicholson-1 ike  scheme  was  developed  and  adapted  to  the 
problem  and  used  to  solve  the  above  set  of  time  and  space  dependent 
differential  equations. 

Plasma  with  parameters  typical  for  a  capillary  discharge  was  found 
to  be  optically  thick  for  resonance  UV  and  for  soft  x-ray  radiation 
at  a  length  over  100  urn;  thus,  to  begin  with,  a  blackbody 
approximation  was  assumed  for  radiation  loss  in  IchKiOO  eV.  A  real 
plasma  radiation  loss  differs  from  a  blackbody  one  (especially  for 
lower  densities)  .  To  correctly  account  for  it  and  to  approach  more 
closely  to  a  real  experiment,  the  factor  for  radiation  loss  was 
introduced  and  is  an  adjustable  parameter  in  the  MHD  calculations. 

The  code  was  tested  to  simulate  Bogen's  experiment  [6).  The 
polyethylene  capillary  they  used  was  2  mm  in  diameter  and  20  mm  in 
length.  The  circuit  capacitor  was  0.3  is F  charged  with  40  kV  maximum 
and  the  circuit  inductance  was  200  nH  (see  Table  1  of  Report  #1  for 
more  details).  The  average  ion  density  was  estimated  as  101*  cm*. 

The  result  of  the  temperature  simulation  is  shown  by  the  solid 
curve  in  Fig.  3.  The  dashed  curve  represents  Bogen's  measurements 
of  temperature  versus  time.  Since  the  experiment  began  from  a 
vacuum  capillary  in  which  the  initial  plasma  density  is  zero,  we  do 
not  expect  our  model  to  simulate  early  development  of  Bogen's 
discharge.  What  we  intend  to  do  here  is  to  simulate  the  cooling 
part.  Figure  3  shows  that  our  rise  time  is  longer  than  in  the 
experiment.  However,  the  result  shows  that  the  plasma  does  cool 
down  from  =  42  eV  to  16  eV  in  about  the  same  amount  of  time.  Thus 
the  experiment  and  the  simulation  have  about  the  same  cooling  rate. 

Figure  4  illustrates  the  power  contribution  of  the  different 
processes.  The  solid  curve  is  E  x  8  heating  power,  while  the  dash 
curve  and  the  dot  dash  curve  represent  the  thermal  cooling  power 
and  radiation  power,  respectively.  From  the  figure,  we  see  that  the 
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electron  heat  conduction  to  the  walls  of  the  capillary  is  about 
triple  the  radiation  cooling.  Milchberg  et  al  [7]  built  a  simple 
model  in  which  they  compared  the  rate  of  thermal  cooling,  radiation 
cooling  and  expansion  cooling  in  a  carbon  capillary  device  and 
found  that  the  thermal  cooling  is  significantly  larger  than  the 
other  two.  Thus  the  thermal  conduction  cooling  mechanism  which 
plays  a  major  role  here  is  consistent  with  Milchberg's  result. 

The  two  total  current  comparisons  show  a  difference  of  less  than  a 
few  percent.  This  is  an  excellent  match  and  it  shows  us  that  the 
magnetic  field  boundary  condition  is  correct  at  relatively  low 
temperatures  (here  less  than  50  eV) .  This  also  gives  us  a  strong 
clue  that  the  code  does  behave  properly  since  the  magnetic  field  at 
the  wall  is  determined  not  only  by  the  external  electric  field  but 
also  by  pressure,  density,  thermal  force,  and  the  degree  of 
ionization . 

The  ion  abundance  evolution  calculated  using  the  general  Saha 
equation  is  illustrated  in  Fig.  5.  In  this  figure,  1  denotes 
neutral  carbon,  2  =  B-like  ion,  3  =  Be-like,  4  =  Li-like  carbon, 
and  5  =  He-like  carbon.  At  the  cooling  time  (from  0.3  to  0.5  tis) 
He-like  ions  dominate  the  abundance  which  is  consistent  with  the 
experiment.  Figure  6  shows  the  plasma  resistance  evolution  curve. 
The  average  value  is  about  1  ohm.  The  ohmic  heating  power  I:R  for 
current  T  =  30  kA  is  about  1  GWatt,  which  is  enough  to  heat  up  the 
plasma  to  45-50  eV  during  a  150-200  nsec  time  interval. 

The  major  advantage  of  our  model  is  that  we  can  monitor  the  1-D 
space  distribution  for  several  quantities.  Figure  7  to  11  show  the 
density,  velocity,  magnetic  field,  pressure  and  temperature 
distribution  as  a  function  of  radius  and  time.  Since  the  electric 
field  E  points  in  the  z  direction  and  the  magnetic  field  B  is  in 
the  6  direction,  the  resulting  force  E  x  B  will  shift  the  plasma 
into  the  center  (pinch-effect).  This  can  be  seen  in  fig.  7.  Figure 
8  shows  that  the  velocity  begins  from  0  and  goes  negative  first, 
which  is  consistent  with  a  z-pinch.  The  fluid  then  sloshes  back  and 
forth,  as  expected  for  MHD  waves.  Because  of  magnetic  confinement, 
and  because  of  the  thermal  conduction  and  radiation  cooling  to  the 
wall,  the  pressure  is  higher  in  the  center  and  lower  at  the  edge  as 
shown  in  Fig.  10.  The  temperature  behaves  the  same  way  but  with  a 
different  slope  due  to  density  gradient,  as  shown  in  Fig.  11. 

The  1-D  MHD  code  has  successfully  passed  several  analytical  tests 
and  was  found  to  correctly  interpret  experimental  data. 

This  code  has  been  used  to  simulate  our  future  experiment  with  Ne~ 
like  Ar.  The  result  of  this  simulation  is  shown  in  Figs.  12-19  for 
2x100  mm  capillary  filled  with  Ar  at  the  initial  pressure  5  mm  Hg. 
Fiq.  12  shows  a  radiated  loss  in  an  Ar  plasma  during  discharge. 
Plasma  resistivity  versus  time  is  shown  in  Fig. 13.  The  circuit 
impedance  VL/C  for  effective  plasma  heating  should  be  choosen 
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close  to  plasma  resistivity,  i.e.  should  be  in  2-5  ohm  range.  Ion 
fluid  velocity  distribution  across  the  capillary  radius  is  shown  in 
Fig.  14.  The  magnetic  field  in  our  capillary  discharge  is  shown  in 
Fig.  15.  Pressure  of  plasma  (which  is  mainly  electron  component 
pressure)  is  shown  in  fig.  16.  Note  essintually  lower  pressure  near 
the  capillary  wall  compared  to  pressure  near  the  center  due  to 
magnetic  confinement  -  this  helps  to  reduce  shock  load  to  the 
capillary  wall  during  discharge.  Fig.  17  shows  Ar  ion  density 
distribution  across  the  capillary.  Ion  abundances  history  versus 
discharge  time  is  shown  in  Fig.  18.  Here  0  represent  neutral  gas, 
8+  -  Ne-like  Ar.  More  than  half  of  total  Ar  ions  are  ionized  to  Ne- 
like  stage.  Temperature  distribution  in  Ar  capillary  discharge  is 
showm  in  Fig.  19.  Temperature  approaches  60  eV  near  the  capillary 
center. 

Currently  the  code  is  under  development  for  simulations  of  a  Ni- 
like  Kr  plasma  capillary  discharge. 


As  shown  above,  a  macroscopic  gain  in  a  capillary  discharge  plasma 
may  be  achieved  for  a  Ne-like  Ar  plasma  with  electron  temperatures 
of  50-60  eV  and  electron  density  of  around  lO19  cm’.  The  calculated 
gain  is  high  for  the  strongest  12-4  transition  (69.76  nm)  and  more 
modest  for  the  three  others.  It  could  be  lower  in  a  real  plasma 
especially  for  plasma  with  non-optimal  parameters.  Hence,  a  long 
length  (10  cm)  capillary  should  be  used.  According  to  the  MHD 
simulation,  plasma  with  the  necessary  initial  density  of  a  neutral 
Ar  may  be  heated  up  to  a  50-60  eV  temperature  and  confined  in  the 
1-mm  radius  capillary  during  several  hundred  nanoseconds  (time 
required  to  get  enough  abundance  of  Ne-like  Ar  at  this  electron 
density) .  The  required  discharge  current  is  a  reasonably  modest  20- 
30  kiloamps.  To  drive  such  a  current  during  hundreds  nanoseconds 
through  1x100  mm  plasma,  a  capacitor  of  0. 1-0.2  (xF  should  be  used 
and  it  should  be  charged  to  40-50  kV.  To  allow  rapid  current  rise, 
the  inductance  of  the  discharge  circuit  should  not  exceed  0.2-0. 5 
M  H. 


IV.  CAPILLARY  DISCHARGE  X-RAY  LASER  DESIGN. 

Using  the  results  of  above  described  simulation  of  Ar  capillary 
discharge,  we  designed  the  experimental  device. 

We  have  began  building  the  capillary  discharge  x-ray  laser  driver 
with  the  following  discharge  circuit  parameters:  C  =  0.22  ft F,  L  = 
0. 1-0.2  nH,  U  =  30-70  kV,  maximum  energy  stored  in  the  capacitor  E 
=  500  J. 


The  capillary  geometry  could  be  varied,  and  for  the  first 
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experiment  a  2  x  100  mm  capillary  was  chosen.  The  length  was  chosen 
longer  than  for  ordinary  capillaries  to  achieve  better 
amplification  for  lines  with  a  lower  gain  coefficient  in  real 
plasma.  The  gas-filled  spark  gap  switch  provides  the  discharge 
triggering.  This  switch  is  attached  to  one  end  of  the  capillary, 
while  another  one  is  open  for  diagnostics.  Two  extra  side  ports  in 
the  capillary  discharge  holder  provide  access  to  the  quarts 
capillary  and  can  be  used  for  observation  of  the  plasma  in  visible 
and  near  UV  light.  The  device  is  very  compact:  6x2x3  dm3  without  the 
capacitor  and  vacuum  pump,  and  even  not  much  larger  with  both.  It 
is  really  table-top  x-ray  laser  driver. 

The  design  of  the  capillary  discharge  is  shown  in  Fig.  20. 


CONCLUSIONS. 

Atomic  and  kinetics  calculations  along  with  MHD  simulation  have 

shown  that  the  proposed  capillary  discharge  plasma  will  serve  as  an 

active  media  for  soft  x-ray  laser  driver.  Calculations  for  Ne-like 

Ar  have  shown  that  a  medium-energy  capillary  discharge  plasma  I 

should  have  parameters  favourable  for  obtaining  ASE  in  Ne-like  Ar, 

with  a  gain  in  the  range  0.3-1  cm'  for  long-wavelength  3p-3s 

transitions  (=70  nm)  and  0.4  cm !  for  shorter  wavelength  (=43  nm) . 

To  explore  such  a  gain,  a  10  cm  length  capillary  with  0 . 5  k J  stored 

energy  is  under  construction  now.  Several  4d-4p  transitions  in  Ni- 

like  Kr  (with  X  =  30  and  40  nm)  have  also  been  found  to  be  the  I 

potential  lasing  lines  for  Kr  filled  capillary  with  lower 

temperature  plasma.  Kinetics  calculation  for  a  gain  estimation  in 

Ni-like  Kr  plasma  and  MHD  simulation  of  a  capillary  discharge 

optimized  for  laser  transitions  in  such  the  plasma  are  on  the  way. 

The  capillary  discharge  x-ray  laser  driver  is  under  construction.  I 


Two  papers  resulted  from  Phase  I  research  efforts  [8,9],  and  more 
expected  to  be  published  during  Phase  II  to  inform  scientific 
community  about  progress  of  our  capillary  discharge  x-ray  laser  • 

project . 
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FIGURE  CAPTIONS. 


Fig.  l.  Calculated  population  of  levels  of  Ne-like  Ar  in  the 

capillary  discharge  plasma  for  Tt=60  eV,  NC=1018  cm3. 

Fig.  2.  Calculated  population  of  levels  of  Ne-like  Ar  in  the 

capillary  discharge  plasma  for  Tc=60  eV,  Nc=1019  cm3. 

Fig.  3.  Temperature  comparison  between  our  simulation  and  Bohen's 
experimental  result. 

Fig.  4.  Power  contribution  of  different  physical  processes  in 

simulation  of  Bogen's  run. 

Fig.  5.  Ion  abundance  in  simulation  of  Bogen's  run. 

Fig.  6.  Plasma  resistance  in  simulation  of  Bogen's  run. 

Fig.  7.  Ion  density  in  simulation  of  Bogen's  run. 

Fig.  8.  Plasma  velocity  in  simulation  of  Bogen's  run. 

Fig.  9.  Magnetic  field  in  simulation  of  Bogen's  run. 

Fig.  10.  Pressure  in  simulation  of  bogen's  run. 

Fig.  11.  Temperature  in  simulation  of  Bogen's  run. 

Fig.  12.  Radiated  power  (W/m3)  for  Ar  plasma  in  simulation  of 

capillary  discharge  x-ray  laser  driver.  Initial  Ar  atomic  density 
is  1017  cm'3,  or  4.8  mm  Hg,  capillary  radius  is  1  mm,  length  100  mm. 

Fig.  13.  Plasma  resistance  (Ohm)  versus  time  in  the  Ar  capillary 
discharge  x-ray  laser  driver. 

Fig.  14.  Ion  velocity  (m/sec)  in  the  Ar  capillary  discharge  x-ray 
laser  driver. 

Fig.  15.  Magnetic  field  (T)  in  the  Ar  capillary  discharge  plasma. 

Fig.  16.  Pressure  (Pa)  in  the  Ar  capillary  discharge  plasma. 

Fig.  17.  Ion  density  (m'3)  in  the  Ar  capillary  discharge  plasma 
simulation. 

Fig.  18.  Abundances  of  different  ionization  species  in  the  Ar 
capillary  discharge  plasma  simulation. 

Fig.  19.  Electron  temperature  (eV)  in  the  Ar  capillary  discharge 
simulation. 


Fig.  20.  Capillary  discharge  x-ray  laser  driver  design. 
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Fig  2.  Level  populations  versus  time  tor  Ne-llke  Ar  In 
a  capillary  discharge  plasma  with  T  ■  60  eV,  N  «  I019cm"3 


abundances 


Fig.  5.  Carbon  ion  abundance  in  simulation  of  Bogen’s  run. 
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Fig.  6.  Plasma  resistance  evolution  curve  in  simulation  of  Bogen’s  run 


0.66 


0.45 


0.60 


Fig-  11.-  Temperature  in  rurnu  Ation  of  Bogen's  run 
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Fig.  20. 


Capillary  discharge  x-ray  laser  design. 


